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I. FUNDAMENTAL PRINCIPLES 
1. Preliminary Remark 



In the older, simple centrifugal impeller theory, it 
is assumed that the flow not only follows the shape of the 
"blades over the blade surface, hut that it approximately 
maintains the same character "between the "blades. This as- 
sumption is justifiable to some extent if the distance be- 
tween the blades is small compared with the radius of cur- 
vature of the absolute streamlines, as is the case for a 
relatively slowly rotating impeller. The more the rota- • 
tional speed was raised, however, and the blade area cor- 
respondingly reduced the more the inadequacy of these 
simple assumptions appeared. Particularly at the blade 
tips, the mean direction of the streamlines no longer 
agreed with the directions of the inlet and exit tangents 
of the blades and these angles had to be corrected by em- 
pirical values. Por very high-speed impellers (Kaplan 
turbines) the simple ideas underlying this method of com- 
putation practically lost their meaning entirely. It was 
then that these blades began to be designed by the pro- 
cedure developed in aeronautics for isolated wings, with 
the effect of the neighboring wings taken into account 
as a disturbance. This method of computation becomes 
simpler the greater the distance between the blades, with 
the resulting small disturbance by the neighboring blades. 
An important aid for this computation procedure is the 
conformal transformation (see reference .1) of a blade casr 
cade. Even by this conformal transformation method, how- 
ever, the computational difficulties increase very greatly 
if the blades are too large or if their distance apart is 
too small. A characteristic that indicates whether a given 
system is suited to the method of conformal transformation 
is that the blades do not overlap (i.e., that it be pos- 
sible to see through between them when viewed along the 
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axial or radial direction for axial and radial impellers, 
respectively). A further disadvantage of the conformal 
transformation method is that it is restricted to two- 
dimensional flows so that it cannot be applied, for exam- 
ple, to impellers of nonparallel casing walls such as that 
shown in figure 1. 

The procedure to be described in what follows is in- 
tended to fill the gaps in the two methods mentioned above. 
We restrict ourselves essentially to radial impellers with 
cylindrical blades since, as Prasil has shown (reference 2), 
the flow about an arbitrarily curved surface of revolution 
(for example, in a Francis wheel) may be reduced to this 
normal form we have chosen by a relatively simple conformal 
transformation. This method starts from the simple hypoth- 
eses of the older centrifugal impeller theory by first as- 
suming an impeller with an infinite number of blades. How 
the flow is modified in passing to a finite number of 
blades is then investigated. The blade shape must be ad- 
justed to this changed flow. For the computation of the 
flow for a finite number of blades, the approximation 
method for isolated wings as developed by Munk, Prandtl 
and Birnbaum, or Glauert (references 3, 4, and 5) is found 
suitable. The essential idea of this method is to replace 
the wing by a vortex sheet and compute the flow as the 
field of these vortices. The shape of the blades is then 
obtained from the condition that the flow must be along 
the surface of the blade. If the blade shape thus obtained 
deviates strongly from that first assumed (on the basis of 
the computation for an infinite number of blades) the com- 
putation is to be repeated with the modified blade shape, 
since the vortex sheets whose field is given by the flow 
should have the same form as the blades. Slight deviations, 
however, are without appreciable effect, as is known from 
the corresponding computations for isolated wings. If the 
blades are not of a negligible thickness, the latter can 
be taken into account in the familiar manner by the intro- 
duction of sources and sinks. 



2. Representation of the Field of Flow 

for Parallel Casing Walls 

It might be possible to proceed by computing the 
field of a vortex which replaces a blade element between 
r and r + dr, adding the disturbance velocities of n, 
such vortices corresponding to the n blades and finally 
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integrating over the "blade length. This computation would 
b*e ! c f ome^ 'rather cumhe ; r.som.e,:.howeYer, since for each point' 
the vortices corresponding to the "blade elements at radius 
r would in general give rise to different disturbance 
velocities. Such computation would become particularly 
tedious when the flow 'cons idered is not two-dimensional 
as "in the case of a wheel with nonparallel walls like that 
of figure 1. The computation simplifies considerably if 
all the point singularities (vortices, sources) are repre- 
sented as vortices of sine wave periodic fluctuation con- 
tinuously distributed over the circle circumference.* Such 
periodic vortex sheets give rise to fields with the same 
periodic fluctuation about the circumference of a circle. 
In figure 2, for example, is shown a periodic vortex dis- 
tribution at radius R and the corresponding distribution 
of the radial and tangential components of the velocity at 
radius r. For characterizing the disturbance only the 
maximum value for each radius need be determined, the de- 
pendence on tho angle cp being directly given by the 
periodic distribution. 

If the velocity difference of the front and back 
sides of the blades at radius E is equal to Av, the 
circulation about a blade element of length ds is equal 
to d T = Av ds. If the element lies between the radii E 
and E + dE, then dE = ds cos d where is the angle 

between the blade tangent and the radius (fig. 3). For n 
blades the n concentrated vortices of circulation d T 
distributed over the circle circumference may be repre- 
sented as continuously distributed vortices, of intensity 
T* per unit length by the Pourier series 

¥ 

T» = j l + 2 Z cos X n (<p - 6) | (l) 

2 E it {_ X=i J 

where 6 = f (r) gives the blade shape. The objection 
might be raised against this representation that this 
series no longer converges and the point vortices can be 
represented only approximately by breaking off the series 
after a few terms. This objection is without justifica- 
tion, however, since in the application of this series to 
the computation of the disturbance ■ velocities the series 
that arise are always convergent. (See equations (7), (8).) 

Tor the idea of replacing the discontinuous blades by 
Fourier series of continuous periodic functions, we are 
Indebted to an oral remark of Mr. J. Ackeret. 
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If the Fourier series (l) is broken off after the 
first term, it would mean that the concentrated vortices 
on the {Wings"Ytare uniformly distributed over the circle 
circumference', i.e., that instead of a finite number of 
blades we assume an infinite number. Therefore, if only 
the first term is taken into account, this would corre- 
spond entirely to the above-mentioned computation proced- 
ure of the older, simple centrifugal impeller theory for 
which the streamline pattern between the blades is assumed 
the same throughout. It is the additional terms of the 
Fourier series that bring out the more accurate conditions 
at the various positions between the blades. The fact 
that in many cases this old, simple theory, that is, the 
single first Fourier term, nevertheless yields useful re- 
sults, leads us to expect that only a very few terms of 
the Fourier series are sufficient to represent the flow 
adequately. The better the agreement of the simple theory 
the fewer the terms that will be needed. A larger number 
of terms therefore will bo mainly required only for the 
neighborhood of the blade tips. In the region of the 
blade tips, however, the neighboring blades play a rela- 
tively small part, so that the conditions at this position 
may be computed on the basis of the knowledge of those 
about the single blade and the effect of the neighboring 
wings taken into account as a correction. 

3. ITonparallel Bounding Walls 

The method of representation of the discontinuous 
vortices by continuous functions shows up to particular 
advantage for nonparallel bounding walls (fig. l). On 
account of the circular symmetry of the bounding walls 
the disturbance field also in this case shows the same 
periodicity characteristics. It is necessary only to com- 
pute the maximum disturbance velocities as a function of 
r/R and n (llumber of periods) for each shape of bounding 
surface. This computation may proceed by first determin- 
ing the relations for the two-dimensional flow. The ef- 
fect of the radially varying distance between the walls is 
to produce an additional increase in the radial velocity 
components v r as a result of the narrowing. The tangen- 
tial components of the velocities v^ are not directly 
affected by the walls since the distance does not vary in 
the tangential direction. They are, however, indirectly 
changed, due to the change in radial velocities by the con- 
dition of freedom from vortices. If h is the distance 
between the walls at radius r and dh/dr, its rate of 
increase with the radius (the inclination of the walls) 
the radial increase in the radial velovity v r at this 
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position is* 

- ■ :■ AJft - (2 ) 

9 r h dr 

This velocity change is superposed on the velocity field 
of the two-dimensional flow in which generally there is 
also a velocity drop 3v rQ /dr. These changes in velocity 

due to the inclination of the hounding walls may also be 
brought out by the introduction of sources and sinks in 
the flow while maintaining the flow two-dimensional. The 
strength of the required additional sources and sinks per 
unit volume is 

av r v r a log h , „ 

E = i- = 1 _ = -v r ■ ■ .. (3) 

9r h dr dr 

The signif icance of these sources and sinks can- also 
be made clear by the following consideration: The two-di- 
mensional flow does not satisfy the boundary condition at 
the walls. The fluid would pass through the bounding walls . 
(fig. 4) with a velocity component normal to the wall v n = 

cL It 

v r sin 8 = -v r cos 6. Through a surface element dP of 

dr 

the bounding wall therefore a quantity Q = v n dP ■ «= <~ 

-v r dP cos 6 passes through, where dP cos 8 is the pro- 
dr 

jection of the surface element in the direction of the axis. 
If vie now apply to this surface element uniformly distrib- 
uted sources of the total strength E = Q, the normal com- 
ponent v n is exactly balanced. There arises only the dif- 
ficulty that the effect of the sources and sinks is not uni- 
formly distributed over the height h if the strength of 
the sources, or v r , in a region which is not large com- 
pared to h, varies appreciably. These sources may not 
therefore be concentrated at the bounding walls but must be 
uniformly distributed over the height h. There is then 
obtained per unit volume the source strength given by equai 
tion (3). 

*We assume here that the distance between the walls is 
small compared with their radius of curvature and that, the 
inclination of the walls is small so that the flow between 
them may be considered constant' on every line parallel to '• 
the axis. If this condition is not sufficiently satisfied 
(for example, in the neighborhood of the axial inlet (fig. 
1)), then the flow may be broken up into layers by the 
method of Prasil and each of the layers separately treated. 
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Since the initial disturbance lying on a circle r 0 
(vortex) of definite frequency gives at every other radius 
a periodic disturbance velocity v r of the same period 
the sources and sinks required for the balancing of the 
nonparallel wall will be distributed according to the same 
period. The added sources and sinks will change the ve- 
locity distribution, but since they are distributed accord- 
ing to the same period as the initial disturbance, only the 
amplitude of the disturbance velocities, but not their 
position and period, will be changed, as already remarked. 

By the above considerations we are in a position to 
compute the flow for an arbitrary shape of the bounding 
walls. Computation of practical examples gives., however, 
the happy result that the shape of the bounding walls in 
almost all practical cases that arise shows up only the 
first term of the Fourier series, that is, in the mean ra- 
dial discharge velocity, whereas, on the higher terms, it 
has only a negligible effect. Since, however, the mean 
radial discharge velocity may be obtained in a quite ele- 
mentary manner from the area of discharge 2 r-rrh, the 
Fourier representation of the source distribution becomes 
superfluous for practical computation. We must, neverthe- 
less, carry out this relatively cumbersome computation in 
order to show that in the majority of cases it is unneces- 
sary. 

II. TEE FIELD FUNCTIONS FOR THE TWO-DIMENSIONAL FLOW 



Following the mothod outlined in the above section, 
we consider the blades to be replaced by vortices and com- 
pute the velocities induced by these vortices, the vor- 
tices dr being periodically distributed according to 
equation (l) on a circle of radius R at the points of 
intersection with the blades. We must now compute the 
velocity field due to such a harmonic vortex distribution. 

Let the vortices lie on the circle of radius R, 
their maximum intensity be p,, the period of distribu- 
tion 2 it/hi, and the maximum values be displaced from the 
zero position by (X + 2-rr/m) (fig. 5). The circulation 
distribution over the circumference is accordingly 

*Y = m- cos m(cp - x) (4) 

where m 1 is an integer. At an arbitrary point (r, cp ) » 
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a tangential velocity v^ and a radial velocity v r are 

produced, the velocities having' the same period as the 
vortex distribution. At the circle of radius B the tan- 
gential components must show a discontinuity, correspond- 
ing to the vortex distribution, of magnitude 



|i cos n(cp - X) 



The radial component v^, from considerations of Continuity, 
can only vary continuously. At infinity and at the origin 
both v^ and v r must tend to zero. Prom these boundary 

conditions and the requirement of a flow free from rotation 
and sources, there are obtained the velocity components 



-m-i 



V <P = 2 CI ) m(0P ~ X) (r > R) 



) cos m(c P ~ (r < B) 

2 * E ' J 



(5) 



and 



v r = 



v r = 



- J±(£) 

2 \R / 
_ H / r 



-m- l 



2 U / 



m- l 



sin m(cp - X) 



sin m(cp - X) 




(6) 



Having determined the velocity field of a harmonic 
vortex distribution, we may now give the velocity field 
due to the discontinuous distribution (1) : 

A Vp (r^p) = - _2L_ L iH Z ( £ ) cos. X n(q> - 6 ) (r< E) 

2itR 9E \ = i \E / 



The part independent of the angle is effective only in the 
external region. Tor the radial velocity, there is further 
obtained 



8 



N.A.C.A. Technical Memorandum Ho. 902 



A v r (r,cp) = - — AE 2 ( L) sinXn(cp - 8) ( r > S ) 

2ttE 3R X=i ^E / 

~.-n oo , _ .. Xn — i 

A v r (r,cp) = - — — AE Z (- J sin Xn(cp - 6) (r<E) 

r ^ 2ttE BE a = i Ve y 

We shall denote the inner blade radius by r^ and the 
outer by r Q (fig. 3). The velocities produced by the 
totality of vortices on the circles between the roots and 
tips of the blades are then given by the relations* 




r 



*In the later computations the velocities always occur mul- 
tiplied by r. 

The series in equations (7) and (8) may be summed. P. 
Staufer (Wass erkraf t u. "tasserwirtsch. 31 (1936) p. 212) 
has also replaced the blades by vortex distributions and 
worked with the finite expressions. This representation 
permits, however, of no generalization to wheels with non- 
parallel bounding walls. 
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III. THE COMPUTATION OF IMPELLERS WITH PARALLEL BOUNDING- WALLS 
WITH THE AID OP THE PERIODIC Pi ELD FUNCTIONS 



The computation procedure is now the following. The 
vortex distribution 3r/dE is prescribed and the corre- 
sponding blade shape is to be determined. There is first 
computed from the previous simple theory the blade shape 
which would correspond to a complete discharge of the water 
or air through the blades, that is, to an infinite number 
of blades. The blades then affect only the tangential ve- 
locities 

r 

n p ar 

rv = — / — dR 
® 2tt J 9R 

' r i 

(this is the first term of equation (?)), i.e., we have 



v = 



~— r(r) 
2ttt 



vrhere T is the circulation of each blade and nT there- 
fore that about all of the n blades. We assume that there 
are no guide vanes, so that the fluid enters without rota- 
tion. In the case of an inflow with rotation of circula- 
tion Tj) the component v rf) would become larger by 

r D /2Trr. Let the mean radial velocity unaffected by the 
blades be c r m » then the relative path is determined by 



(9) 



(fig. 3) where # denotes the angle of the relative flow 
path with the radius (for backward curved blades $ is 
always negative (fig. 3)); w is the angular velocity of 
the impeller. Since it is convenient to compute with non- 
dimensional magnitudes, we introduce 



T*(r) 



T(r) 



2-rrr a u> r a 



v* = 
v r 



id r. 



v* = 

CP 



(JU 



L 
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Equation (9) then "becomes 




(9a) 



from which the mean flow path may readily he computed and 
found to he 



so that 



tan d 



cp (r) - cp(ri) = 



rdcp 
dr 



t an $ ■ 
a r 



For an infinite number of blades the mean flow path and 
blade shape agree so that the function 6(r) (fig. 3), 
which gives the blade shape, is also determined. 

With the aid of equations (7) and (8) it is now pos- 
sible to compute the additional velocities for a finite 
number of blades. Since, in the cases that practically 
arise term by term, integration is possible, we compute by 
the modified formulas 



-H-v* = nF*(r)+n £ 

r a cp 



r 

co a 



ar /„ , - An 

cos Anfo, *■ 6 )- d E 



- e ) da 



(7a) 
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r a r 



: • ■ CO. 

•n £ 



(I) 



sin Xn(cp - 6) d r 



(8a) 



sin Xn(cp - 6) d R 



and 



— v* = — v* - nr *(r) 

r a cp add r a co 



r r add 



r v * 



Assuming that 9r*/3 r for r^ and r a , that is, "blade 

root and tip, tends toward zero, the convergence is 
assured and rapid. On the proper choice of 9 P/Br , we 
shall have something. to say later. 

After v^ and v* a ^ have teen computed, the 

new blade shape is determined 

6 x (r) - 6 l(ri ) = J IfiiJjL dr 



with 



"•^ + add " (^) 



tan = (10) 



c*_ ( - 



This deviates from the old shape, particularly at the root 
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and tip. The velocity field is therefore computed anew 
"by assuming the concentrated vortices on the new blade 
shape, i.e., in equations (7a) and (8a) replacing 6 "by 
b x . This method is continued until 6 n agrees with 

6 n _ 1 . Practically & 3 will already deviate little from 

With regard to the numerical evaluation, it is to be 
noted that, on account of the high exponents Xn, the in- 
tegrands occurring in (7) or (8) have appreciable values 
only about the point considered (r,cp), so that in the 
greatest portion of the blade region only few Fourier terms 
will be required. Only at the blade tips where the vortex 
distribution drops very rapidly is it necessary to take a 
relatively large number of terms of the series (7) and (8) 
in order to obtain accurate values of the additional veloc- 
ities at the blade root and tip. This difficulty may.be 
evaded in the following manner: The integrals are evalu- 
ated, for example, if we consider the blade tip only, up 
to a value r, < r 0 where r, is so chosen that the 
blade between r x and r a is as nearly as possible re- 
placeable by a straight line. For the vortex distribution 
between r, and r a , there is determined the center of 
gravity, which we assume is at r^. The tips of all blades 
except the one under consideration are now replaced by a 
single vortex and the effect of the remaining blades is 
thus computed. At the blade itself th e circ ulation distri- 
bution is replaced by the function c vr a - r as may cer- 
tainly be done with good approximation and thus the addi- 
tional velocities are computed. The constant c must be 
so determined that a good continuity is obtained with the 
circulation distribution in the region r< t 1 . This pro- 
cedure has also been applied to the example considered be- 
low for shortening the computation. 



IV. IMPELLERS WITH NONPARALLEL BOUNDING WALLS 
1. The Impeller Width Decreases According to a Power Law 



In part I, 3 it was explained that for the case of 
nonparallsl walls there are additional sources whose 
strength is given by the relation 



v r dh 
h dr 



(3) 
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In the above equation r r is composed of a portion v r 

1 

corresponding to the discharge flow, a portion , v_ cor- 

1 2 

responding to the vortex flow, and a portion v r due 

to the sources themselves. The determination of the 
sources leads, for known v™ and v- . to the solution 

of an integral equation. For the case that the width var- 
ies by a power law, that is, where the width is given by 

h = C r~S (11) 

the computation of the flow may be more simply carried 
out. We may avoid the solution of the integral equation 
if we proceed as follows. 

We seek to determine a periodic velocity field which 
everywhere has sources, that is, whose divergence is dif- 
ferent from zero, and which, in addition, possesses on a 
circle of prescribed radius a discontinuity in the tangen- 
tial velocity corresponding to a given periodic vortex 
distribution. By superposition then may be determined the 
velocity field which has vortex distribution of different 
periods on many circles, as corresponds to the circulation 
distribution of the impeller under consideration. 

Lot the unknov/n radial and tangential velocities have 
the period m and amplitudes v rQ and v^ 

v r = v r 0 sin mrp v cp = T cp cos m CP (12) 

(See equations (5) and (6).) lor the field divergence 
which is then also a periodic function, there is obtained 

div v = E m sin mco = j^r( rv r 0 sin mcp)+~(v r p 0 cos m <$ J (13) 



Carrying out the differentiation and dividing by 
we obtain for the amplitude of the source 

E m = ij A ( r v r ) - m v 1 
m r L Sr r o cpo J 

Since the required flow is free from vortices, except on 
the single circle of radius E that has a prescribed vor- 
tex distribution, we have for r < E and r > E 



sin m cp 

(14) 
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rotv=0=— — — (r v C n cos mro) - ^ (v r „ sin m cp) 
— r or 0 rSrp ° 

or 

v„ + r — m v r = 0 (15) 

3r 0 

From equations (14) and (15) there is obtained by elimina- 
tion of the tangential velocity and its derivatives 

2 Effi + r ^III = (1 _ »») Ho + 3 ^£o + r ^ZEJL (16) 
3r r Br Br 2 

In order that the sources correspond to the prescribed de- 
crease in the width of the impeller, they must satisfy the 
condition (see equation (3)) 

, _ V *" 0 dh 

Em — r 17 

For the case where the width h decreases to a negative 
power h = C r~&, there is obtained 

E m = f v r 0 <") 



We have then from equation (16)* 



2 



tr / * 9 T r e \ v r a v r 3 v r 

r r o V r V a r r s r °y r 3 r 9r s 

or after rearrangement 



2 

L_lEa + Llla (3 - g) + ^-°-[(l - ma) - g] = 0 (18) 
3 r 3 3 r r 



*It may be readily seen that every other function h' 
makes the solution of the differential equation difficult. 
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This differential equation is solvable "by an exponential 
substitution . .. . 

v = a r K 
r o 

where a and 0 are constants still to he determined. 
For the constant P, there is obtained 



i.e., there is a solution for the amplitude of the radial 
velocity 



(vr o )i = ai r ' = a ± r Pi (19) 



which vanishes for r — > 0, and a solution 



- C 1 - §y</^ 



£1 

4 



which is zero for r > oc Both of these solutions must 

then be so assembled that the corresponding velocity field 
has exactly the prescribed discontinuity in the tangential 
velocity field on the circle with radius R. Prom equation 
(14), there is obtained for the amplitude of the tangential 
velocity field 



° V % = v r 0 + * - T ~ - r E m 

and, further., making use of equation (17) 

d v r 

m v o = ^rA 1 - g) + r — - — - 
^o A o 3 x 

For the difference of the amplitudes of the tangential 
velocities on the circle with radius R, we obtain after 
substituting the solutions (19) and (20) 
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- V 



o a 



cPo i m 



= i (a a 0 a H Pa - a ± Pi 



Since at the position of the discontinuity of the tangen- 
tial velocity, we have on account of the continuity of 
the normal velocity 

a a R Pa = ai E Pi 



there is obtained , 

, IT7 -(i-^V /«.=+ 4- 



C0 o a °Po i m 
The difference v 



m 



- V 



o a 



-a a E 



represents the amplitude of 



o 1 



the vortex distribution 



- v 



2 n / ar 



^o a r 'o i 2ttE v 9E 
where = m. There is then obtained the constant 



a. 



2 n / ar 



nX 



E 



2ttE n 3E ' /" 22 , ^ 
V 4n" x + g 2 



- ^ 4E ) 1 

V SR / / 

7 1 



2ttR V 3E 



E 



'1 + 



4 n 2 X 3 



(l - |) + nx/T7 



?> (2 



4 n 3 X 2 



III! III! Ml lllllllllim Mil ■ II I III III! I IH ■■■■■■■ 
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Jor the velocities in the radial direction, we obtain 



, _ _L./9r..\ i 

r a~"i^sVaH AR ;~7= 



el 

4 



sin Xncp 



31 



4n 2 X 2 



1 2ttR J ) r~ v S.J 



2 + li 



>(22) 



sin Xnco 



4n 3 X 2 

and for the tangential direction 



2ttR VAB / 



a. 2ttE ^3R 



cos Xnci 



4n2Xs 



> (23) 



2ttR v 3R ' 



2 V A 4/r\ x 

^— J cos Xncp 



e 



— V R' 



4n 2 V 



There was thus obtained, in place of (7) and (8), for an 
impeller width decreasing according to a power law, for a 
finite number of blades 
r 



2tt 



9£ 
dR 



1+ Z 2 - 



4 -(~) cos Xn(q>-e) 



d R 



► (24) 



a 

_n_ r-_ar 

2tt / 3R 
J r 



Z 
X=l 



«. JT^S. Pi+i 



4 /r 



(jgj cos >n((p-6) 



d R 
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and the additional radial velocity 
r 



n / aF v 
r v r = - — / — £ 



2tt 



SR >,=i 



\Pa +1 

1 (^Lj sin Xn(q)-6) d E 



1 + 



e 2 



4 n 2 X 3 



= (i) 




4 n 2 X« 



sin X n(cp - 6) d E 

(25) 



These relations (24) and (25) for hyperbolically decreasing 
impeller width are employed in place of (7) and (8). Other- 
v/ise, the computation proceeds in the same manner. The 



mean radial velocity c- 



m 



for decreasing width is naturally 



different from that for a constant width, 



It is to be noted that 




nV + 



g 2 



(26) 



and, hence, for large nX 



Pa + 



1 = - nA 



3 i + l = + n> > 



£ _ 



S 



• • . i 



(27) 



From this approximation, it follows that the effect of the 
decreasing width on the harmonics for impellers of "blade 
number > 2 is in general very small, so that the variable 
width need be taken into account only for mean radial dis- 
charge velocities c r = Q/2rTrh, but otherwise the compu- 

m 

tation may proceed as though the walls were parallel. 



N.A.C.A, Technical Memorandum No. 902 



19 



2. Arbitrary Decrease in Impeller Width 

For an arbitrary decrease in the impeller width h(r) 
there are obtained in place of (18) differential equations 
that no longer have a simple solution^ We shall therefore 
determine indirectly the sources required for satisfying 
the boundary conditions. The source strength B is given 
by the equation 

_ v r dh 1 dh . v # — „\ 

E = - = (v r + v r + v r ) (28) 

h dr h dr r i r z r » 

where v r corresponds to the pure discharge velocity: 

v_ = c_ . v_ is induced by the vortex flow and is ob- 
1 l r m r s 

tained from equation (8) so that it is a periodic function 
in cp and the portion v r3 is the radial velocity due to 

the sources themselves. To solve equation (28) there is 

assumed for the unknown sources a periodic expression of 
the form 

r)+2 Z a^ n (r) cos Xncp + *>Xn( r ) sin * n cpj (29) 

with the undetermined coefficients a, (r) and b. (r). 

An An 

The radial velocity induced by these sources is similarly 
a periodic function in cp J 
r r 

- An , -An 



n 

V r, = 3irr "* 




J a 0 (R)+ ^Ja^ n (E) (|) cos X ncD+b An (R) (§) 



sin An cp 



, ■ An 
A n 



d a 



(R) (5 ) cos An(p+b^ n (R) ^ sinXncri 



dS 



(30) 



Substituting the values v r , v r , E, and v r in equation 

(28), the latter breaks up into an infinite number of equa- 
tions corresponding to the number of harmonics and there is 
obtained an infinite set of independent integral equations 
for the unknown coefficient functions a An ( r ) &n & ^An^ r ^* 

These integral equations may be solved by iteration. We 
shall not, for the reasons given above, go into the solution 
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of these integral equations, since it is only very rarely 
that it becomes necessary to solve them. 



V. COMPUTATION OF AH IMPELLER 



In order to illustrate the method of computations of 
centrifugal impellers, taking into account the finite num- 
ber of "blades and constant or decreasing impeller width, 
we shall now consider an example. Let it be required to 
design a centrifugal impeller (blower) which imparts to a 
discharge of 600 m 3 per minute a pressure rise of 200 mm 
water (air conditioner for mine operation), that is, 



Q = 600 m 3 /min = 10 m 3 /s ^ 

) (31) 



P tot = 200 mm WS = 200 kg / m 

Let there be available a motor drive with 720 r.p.m., cor- 

2-nn _ , 

responding to an angular velocity u) = qq = 75.5 s . The 

suction velocity may be assumed as c s = 25 m/s. If the 
diameter of the impeller at the inlet is made equal to the 
diameter of the suction pipe, there is obtained from 

r? tt c s = Q for the inlet radius r^ = 0.36 m. If it is 
assumed in addition that the suction and inlet velocity are 
approximately the same ( c r m . = c s = 25 m/s) there is ob- 
tained for the width of the impeller at the inlet from the 
relation c r m . x 2 r^ tt h^ = Q the value h^ = 0.18 m. 

The outer diameter of the impeller is chosen as 2 r a = 

1.4 m, that is, about equal to twice the value of the in- 
ner diameter. For constant impeller width the mean radial 
exit velocity then becomes 



c r = 12.84 m/s(c£ = c„ /iur. = 0.243) 

The desired increase in pressure determines the required 
total circulation* 



*In relation (32) we have assumed that the efficiency n = 1 
so that friction losses are not taken into account. If an 
efficiency of 75 percent is assumed, the impeller would de- 
liver a pressure head of 150 mm of water. 
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A p = p m £*L (32) 

.'..,-..„.., 2 71 



fp = 1 kg sec 3 \ ^ projn a p = 200 kg/m 2 , there follows 
\ 8 m ' 
for the nondimensional total circulation 

n r a td . 

n T* = = — = 0.572 (33) 

2 7t r a ou r a p r| 

With these values, there is obtained from (9) for the case 
of an infinite num"ber of "blades, for the "blade tangent at 
the inlet 

tan * 4 = = 1 = -1.1 (34) 

Cr *i Q 

2 ti hi 

that is, the angle $ ^ = - 47.5°, and, from (9) and (9a), 
for the blade tangent at the exit for constant wheel width 

n r 

- <d r. 



2 it r a a n T* - 1 , . /„,->. 

tan - a = _ = - 1,8 (35) 

c r c* 
n a . r m a 



i.e., 



4 = -60.5° 



The absolute inlet velocity c ri amounts to 25 m/ s (see 
above) and the relative velocity 



c r . 



v rel , =. " X1 = • = 37 m/s (36) 

1 cos ^ COS ^ 

The relative exit velocity is obtained as 

c r . 

v rel = = 26 m/s (37) 



a cos J a 



and the absolute* 



*It should be noted that for backward curved blades £ is 
negative. 
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0 exit 



= / v rel n + 0) s rj + 2 v rel u> r & sin * a (38) 



or' 



3* 



We shall design the wheel with 10 blades. In order 
to determine the blade shape, we must first prescribe the 
variation of the circulation along each blade. Since we 
are considering impact-free entry, the circulation at the 
root and tip of the blade must drop to zero. Along the 
length of the blade the circulation may be assgned at will. 
Having assigned the circulation along the blade, we deter- 
mine the pressure variation, \irhich is of decided importance 
for the behavior of the boundary layer and its separation. 
Since, at the present time, however, these questions are 
not yet sufficiently cleared up, we shall be satisfied with 
an estimate of the blade loading on the basis of a compari- 
son with airfoils. In order to obtain approximately uni- 
form loading of the blades, we have assumed the distribution 
n d r*/ar, shown in figure 6. The distribution of the 
circulation itself is then given in figure 7. With the aid 
of equation (9a) and the one obtained from it there were 
then drawn the relative flow paths for the case of an in- 
finite number of blades (fig. 8). Tor an assumption of 10 
blades, there would thus be obtained as a rough approxima- 
tion (computed as though for an infinite number of blades) 
the blade passage shown in figure 8. The circulation 
about each blade amounts to 0.0572 x 2 tt r a x u; r a . If 

the lift coefficient of the blade is estimated on the basis 
of airfoil investigations, we have 

e = i_£ = 4 tt^-^T* = 4^^^ 0.0572 (40) 
a v t t V t V 



It is not clear which value for 7 must be substituted 
since the formula is itself applicable only to a wing in 
parallel flow. If we assume as a basis the mean value 
between the relative inlet and outlet velocity, we obtain 
c ~ 1.57. If, however, the value at the flow-off edge 
is considered as the determining velocity, there is ob- 



*This equation can be written down immediately or obtained 
from (38) with the aid of (35). 
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tained c a •» 1.9. The e a value is rather high but the 
chbsoh"di'mensions of the--wh«el still may he considered as 
suitable. . 

The radial and tangential additional velocities that 
must he taken into account for a finite number of blades 
were determined with the aid of equations (7a) and (8a) and 
then from equation (10) the new improved blade shape S x 
was computed (fig. 9). Since the shape S x deviated con- 
siderably from the shape S 0 , there was determined a still 
further improvement S 2 (fig. 9). This shape lies so 
close to S x that a further continuation of the iteration 
was dispensed with. 

The pressure distribution at the blade, a factor of 
importance for the efficiency of the impeller insofar as 
the nonaccounted-f or effects of the boundary layer may be 
judged from it, can be determined with the aid of the 
Bernoulli equation extended to rotational flows 

P + f v rel = const + f u,2r3 (41) 

The constant is P^ot- ~ |r c s + Po» if we denote the 

static pressure at the entry to the wheel by p Q . The 
relative velocities at the pressure and suction sides of 
the blade, which is assumed as very thin, differ at the 
corresponding points by the amount of the vortex distri- 
bution there, d^/ds (where s, the arc length of the . 
blade, is measured from the inlet edge). Denoting the 
mean value of the relative velocities at the pressure and 

suction sides by v„ , , we obtain from equation (41) 

rei ra 

for the pressure at the suction side p g g 



P sg 



pt. tl -f-v-f + .(«) 



and for that at the pressure side p dr 

Figure 10 shows the pressure distribution at the impeller. 
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Since, in general, impellers are designed with de- 
creasing width in order to obtain as constant a mean ra- 
dial velocity as possible, we shall also compute the ex- 
ample for the same discharge quantity and pressure rise 
for an impeller with decreasing width. 

Let the width decrease linearly with the radius from 
inlet to exit, so that the exit width h a = 0.11 m (fig. 
11). From the continuity relation 

v r a 2 * r a h a = v r ± 2 " r i h i 

there is then obtained for the mean radial exit velocity 

v r = 0.84 X 25 = 21 m/s 

a 

First we compute the flow paths for the case of infinite 
number of blades. The computation proceeds exactly as for 
parallel bounding walls except that the change in the mean 
radial discharge velocity as a result of the decreasing 
width must be taken into account. We obtain the mean flow 
path indicated in figure 12 by S a . We now take account 
of the finite number of blades by computing the harmonic 
additional velocities from equations (7a) and (8a) where 
we may leave the decreasing impeller width out of account. 
By carrying along two terms of the series, we obtain the 
blade shape S a . To judge the suitability of the blade 
shape, we determine the load coefficient from equation (40). 
We obtain the value 2.1 if we employ the exit velocity 
as a basis for the approximation. The value has thus not 
dropped in spite of the increase in the mean discharge 
velocity. This is due on the one hand to the smaller back 
flow of the blades and on the other hand to the shortening 
of the chords. 

In order to find the effect of the decreasing impeller 
width on the harmonic additional velocities according to 
the discussion of section IV, we computed the additional 
velocities for the case of equal total decrease of width 
but with the latter varying hyperbolically (h = C r~S = 
0.085 r~ °* 74 ) . The difference of the correction values 
for v r and v^ in the harmonics was very small, as com- 
pared with the linearly decreasing width, where no account 
was taken of the decreasing width in computing the harmon- 
ics. Thus there was obtained for the integral value of the 
first harmonic in the formula for the circumferential ve- 
locity (equation (24) made nondimensional corresponding to 
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(7a)) the value 0.0332 compared to 0.0337, the value that 
was obtained without taking account of the decreasing width 
for the harmonic additional velocities. Since the fourth 
decimal, is not quite certain, the value is practically the 
same. Similar relations were obtained for the radial ve- 
locities. Since the effect of the decreasing impeller 
width must he large for a small number of blades, we also 
treated the problem for a wheel with 5 blades and one with 
2. In the table below are given the integral values of the 
first harmonics of the circumferential velocity. 



No. of 


First cog ine 


correction 


First sine correction 


blades 


Decreasing 
width 


Constant 
width 


Decreas ing 
width 


Constant 
width 


10 


0.0332 


0.0337 


0.053 


0.050 


5 


,1009 


. 1004 


.119 


.110 


2 


.306 


.286 


. 164 


. 145 



It may be seen that only for an impeller with 2 
blades does any difference occur that is greater than the 
computational accuracy. It nay therefore be said that it 
is not necessary in the case of impellers the number of 
whose blades is considerably larger than 2 to take account 
of the variable width in the computation of the harmonic 
additional velocities. This may also have been expected 
from a consideration of equations (26) and (27). 



VI. SUMMARY 



For centrifugal impellers with constant and with de- 
creasing width, a procedure was developed which makes pos- 
sible the computation of the unknown blade shape for a 
prescribed circulation distribution over the circumference. 
It was found that the decrease in the impeller width es- 
sentially affects only the mean radial discharge velocity. 
The procedure fills the gap between the two often consid- 
ered limiting cases, namely, that in which the impeller 
is computed as though for an infinite number of blades 
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(#Uler equation) and that in which conformal transforma- 
tion is applied* a procedure that is convenient only in 
the case of a small number of blades. Since as yet only 
the case of impact-free entry has "been treated, only the 
theoretically test operating condition can he computed. 
As soon as, for example, the discharge quantity varies 
the entry angle no longer agrees and every "blade is sub- 
ject to the flow at the entry edge. The loss that is 
thereby incurred may be estimated with the aid of the 
methods of the airfoil theory. This will be discussed in 
a succeeding paper. 



Translation by S. Reiss, 
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Pigs. 1,2,3,4,5,6,7 



Figure 1.- Impeller with 
variable width. 





Figure 2.- Periodic vortex distri- 
bution on the circle of 
radius R and corresponding tangential 
and radial velocities v ? and v r on 
the circle of radius r. 




Figure 3.- Notations. 



Figure 4.- Notations for impeller 
with variable width. 




Figure 5.- Harmonic vortex distribution 
on circle of radius R. 




Figure 6,- Circulation distribu- Figure 7.- Circulation as a 
tion in radial function of the 

direction. radius. 
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Figs. 8,9,10,11,12 




& a =-60.5° 
£ "=-53.0° 
1 =-46.30 




Figure 8.— Blade shape for the assumption 

of infinitely many blades Figure 9.- Taking account of 

' ' (constant impeller width). " the finite number 

of blades in the computation 
of the blade shape. Sq initial 
shape (see fig. 8), S]_ and Sg 
first and second improved 
shapes respectively (constant 
blade width). 
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Figure 11.- Decrease of the 
impeller width. 
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Figure 10.- Pressure distribution 

at the blade Sg 
The pressure which 
corresponds to the 
mean velocity 
is shown dotted. t 




a = -47.6° 
ai " - 39 ' 5 ° 



Figure 12.- Blade shape for equal circulation distribution (fig. 6) 

but decreasing impeller width. S a shape according to 
the theory for an infinite number of blades, S ai shape taking into 
account the finite number of blades. 
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